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Abstract. We present a theoretical model that explains the high energy phenomenology of the 
neighborhood of SNR IC 443, as observed with the Major Atmospheric Gamma Imaging Cherenkov 
(MAGIC) telescope and the Energetic Gamma-Ray Experiment Telescope (EGRET). We also 
discuss how the model can be tested with observations by the Fermi Gamma-ray Large Area Space 
Telescope. We interpret MAGIC J0616+225 as delayed TeV emission of cosmic-rays diffusing from 
IC 443 and interacting with a known cloud located at a distance of about 20 pc in the foreground 
of the remnant. This scenario naturally explains the displacement between EGRET and MAGIC 
sources, their fluxes, and their spectra. Finally, we predict how this context can be observed by 
Fermi. 
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INTRODUCTION 

Recently, MAGIC presented the results of observations towards SNR IC 443, yielding 
to the detection of a new source of y-rays, J06 16+225 (Albert et al. 2007). This source 
is located at (RA,DEC)=(06 h 16 m 43 s , +22°31'48"), with a statistical positional error of 
1 .5' , and a systematic error of 1 ' . A simple power law was fitted to the measured spectral 
points: dN y /(dAdtdE) = (1.0 ±0.2) x 10~ n (£/0.4TeV) _3 - 1±a3 cnrV^eV -1 , with 
quoted errors being statistical. No variability was found along the observation time (over 
one year). No significant tails nor extended structure was found at the MAGIC angular 
resolution. 

MAGIC J06 16+225 is displaced with respect to the position of the non- variable 
(Torres et al. 2001) EGRET source 3EG J0617+2238 (Hartman et al. 1999). Indeed, 
the EGRET central position is located directly towards the SNR, whereas the MAGIC 
source is south of it, close to the 95% CL contour of the EGRET detection. As Albert 
et al. (2007) showed, the MAGIC source is located at the position of a giant cloud in 
front of the SNR, it would not be surprising if they are related, which we explore here. 
The EGRET flux is (51.4±3.5) xl0~ 8 ph cm~ 2 s , with a photon spectral index of 
2.01 ±0.06. Extrapolating the spectrum of the EGRET source into the VHE regime, 
we would obtain a higher flux and a harder spectrum than which was observed for 
MAGIC J06 16+225, supporting the view that a direct extrapolation of this and other 
EGRET measurements into the VHE range is not valid (Funk et al. 2008). 

Here we present a theoretical model (see Aharonian & Atoyan 1996, also Gabici & 



Aharonian 2007) explaining the high energy phenomenology of IC 443, making focus in 
the displacement between EGRET and MAGIC sources. Our interpretation of MAGIC 
J06 16+225 is that it is delayed TeV emission of cosmic -rays (CRs) diffusing from the 
SNR. Finally, we discuss how the model can be tested with observations with the Fermi 
Gamma-ray Large Area Telescope. 

THE SNR IC 443 IN CONTEXT OF THE MODEL 

IC 443 is one of the most studied SNR region at all frequencies. It is an asymmetric 
shell-type SNR with a diameter of ~45 arc minutes (e.g., Fesen & Kirshner 1980). Two 
half shells appear in optical and radio images (e.g, Braun & Strom 1986, Leahy 2004, 
Lasker et al. 1990). The interaction region, with evidence for multiple dense clumps, is 
also seen in 2MASS images (e.g. Rho et al. 2001). In radio, IC 443 has a spectral index of 
0.36, and a flux density of 160 Jy at 1 GHz (Green 2004). Claussen et al. (1997) reported 
the presence of maser emission at 1720 MHz at (/, b) ~ ( — 171.0, 2.9). Recently, Hewitt 
et al. (2006) confirmed Claussen's et al. measurements and discovered weaker maser 
sources in the region of interaction. IC 443 is a prominent X-ray source, observed 
with Rosat (Asaoka & Aschenbach 1994), ASCA (Keohane 1997), XMM (Bocchino 
& Bykov 2000, 2001, 2003, Bykov et al. 2005, Troja et al. 2006), and Chandra (Olbert 
et al. 2001 , Gaensler et al. 2006). The works by Troja et al (2006) and Bykov et al. (2008) 
summarize these observations. In what follows we present some additional features of 
IC 443, relevant for our model. 

Age: IC 443 is agreed to have a middle-age of about 3 x 10 4 yrs. This age has been 
initially advocated by Lozinkaya (1981) and was later consistently obtained as a result 
of the SNR evolution model (Chevalier 1999). 

Distance: Kinematical distances from optical systemic velocities span from 0.7 to 1.5 
kpc (e.g., Lozinskaya 1981). The assumption that the SNR is associated with a nearby 
HII region, S249, implies a distance of ~ 1.5 — 2.0 kpc. Several authors claimed that the 
photometric distance is more reliable (e.g, Rosado et al. 2007), and concurrently with 
all other works on IC 443, we adopt its distance as 1.5 kpc (thus, 1 arcmin corresponds 
to 0.44 pc). 

Energy of the explosion: There is no clear indicator for £51, the energy of the explo- 
sion in units of 10 51 erg. Lacking a strong reason for other numerical assumption, we 
will assume that £51 = 1, although to be conservative, we will subsequently assume that 
only 5% of this energy is converted into relativistic CRs. Reasonable differences in our 
assumed value of £51 are not expected to have any impact on this model. 

The molecular environment: Cornett et al. (1977) and DeNoyer (1979) were among 
the first to present detailed observations of molecular lines towards IC 443. Subse- 
quently, Dickman et al. (1992), Seta et al. (1998), Butt et al. (2003) and Torres et al. 
(2003) among others, presented further analysis. These works conform the current pic- 
ture for the environment of IC 443: a total mass of ~ 1.1 x 10 4 M mainly located in 
a quiescent cloud in front of the remnant (with linear scales of a few parsecs and den- 
sities of a few hundred particles cm~ 3 ) that is absorbing optical and X-ray radiation 
(e.g., Lasker 1990, Troja et al. 2006), a scenario already put forward by Cornett et al. 
(1977). Dickman et al. (1992) estimated that 500-2000 M are directly perturbed by 



the shock in the northern region of interaction, near the SNR itself. Huang et al. (1986) 
found several clumps of molecular material along this interacting shell, with subparsec 
linear scales. Rosado et al. (2007) found inhomogeneities down to 0.007 pc. As it is 
usual, we will neglect these latter inhomogeneities when considering the propagation of 
CRs in the ISM, i.e. we thus assume an homogeneous medium of typical ISM density 
where CRs diffuse. Then, the molecular mass scenario is a main giant cloud in front of 
the SNR containing most of the quiescent molecular material found in the region, and 
smaller cloud(s) totalizing the remaining mass located closer to the SNR. 

DIFFUSION OF CRS FROM IC 443 

The spectrum of y-rays generated through 7T°-decay at a source of proton density n p 
is Fy(E y ) = 2 JJi,, (F^Ex)/ \/El-ml) dE K , where E™ m (E y ) = £ y + m|/4£ y , and 

F„(E„) = Ann p j ' ^ n J p (E)(do n (E n , E p )/dE n ) dE p . Here, da K (E n , E p )/dE K is the 

differential cross-section for the production of 7T°-mesons of energy E K by a proton 
of energy E p in a pp collision. For an study of different parameterizations of this cross 
section see Domingo-Santamaria & Torres (2005) and Kelner et al. (2006). The limits 
of integration in the last expression are obtained by kinematic considerations (see e.g., 
Torres 2004). In these expressions we have implicitly neglected any possible gradient of 
cosmic-ray density in the cloud as well as in the cloud's gas number density. 

The CR spectrum is given by J P (E, r, t) = [cj8/4^]/, where f(E, r, t) is the distribu- 
tion function of protons at an instant t and distance r from the source. The distribution 
function satisfies the radial-temporal-energy dependent diffusion equation (Ginzburg & 
Syrovatskii 1964): (df/dt) = (D(E)/r 2 )(d/dr)r 2 (df/dr) + {3 /BE) (Pf)+Q, where 
P = —dE/dt is the energy loss rate of the particles, Q = Q{E, r, t) is the source 
function, and D(E) is the diffusion coefficient, for which we assume here that it de- 
pends only on the particle's energy. The energy loss rate are due to ionization and 
nuclear interactions, with the latter dominating over the former for energies larger 
than 1 GeV. The nuclear loss rate is P nuc = E/x pp , with x pp = {n p cKC pp )~ l be- 
ing the timescale for the corresponding nuclear loss, k ~ 0.45 being the inelastic- 
ity of the interaction, and o pp being the cross section (Gaisser 1990). Aharonian & 
Atoyan (1996) presented a solution for the diffusion equation for an arbitrary energy 
loss term, diffusion coefficient, and impulsive injection spectrum f m j(E), such that 

Q(E,r,t) = Nof m j(E)8r 8 (t). For the particular case in which D(E) <x E s and/i n j °cE~ a , 

the general solution is f(E, r,t) ~ (N E- a /n 3 / 2 R 3 dif ) exp [-(a - l)t/x pp - (R/R iif ) 2 ] , 

where = 2(D(E)t [exp(/5/ x pp ) — 1] / [tS / x pp ] ) 1//2 stands for the radius of the sphere 
up to which the particles of energy E have time to propagate after their injection. In case 
of continuous injection of accelerated particles, given by Q(E, t) = QoE~ a 2F{t), the 
previous solution needs to be convolved with the function 3T(t — t') in the time interval 
< t' < t. If the source is described by a Heavside function, ST{t) = ®(t) Atoyan et al. 
(1995) have found a general solution for the diffusion equation with arbitrary injection 
spectrum, which with the listed assumptions and for times t less than the energy loss 
time, leads to: f(E, r, t) = {QqE^ /AKD{E)r){2/ ^/k) H r e'^dx. We will assume 



that a = 2.2 and make use of these solutions in what follows. 

Fig. 1 shows the current CR spectrum generated by IC 443 at two different distances 
from the accelerator, 10 (solid) and 30 (dashed) pc. The SNR is considered both as a 
continuous accelerator with a relativistic proton power of L p = 5 x 10 37 erg s _1 (the 
proton luminosity is such that the energy injected into relativistic CRs through the SNR 
age is 5 x 10 49 erg), and an impulsive injector with the same total power (injection of 
high energy particles occur in a much shorter time than the SNR age). The horizontal 
line in Fig. 1 marks the CR spectrum near Earth, so that the excess of CRs in the SNR 
environment can be seen. For this example, the diffusion coefficient at 10 GeV, D\q, 
was chosen as 10 26 cm 2 s , with 8 = 0.5. CRs propagate through the ISM, assumed 
to have a typical density. In the scale of Fig. 1, curves for hjsm = 0.5, 1,5, and 10 cm -3 
would be superimposed, so that hjsm becomes an irrelevant parameter in this range (this 
stems from the fact that the timescale for nuclear loss % pp obtained with the densities 
considered for the interestellar medium, hjsm, is orders of magnitude larger than the age 
of the accelerator). Differences between the different kind of accelerators assumed are 
also minimal for the SNR parameters. 

Fig. 2 shows the result for the y-ray emission coming from the cloud located at 
the position of the MAGIC source, when we assume it lies at different distances in 
front of IC 443. The giant cloud mass is assumed (consistently with observations) as 
8000 M . The accelerator properties and power of IC 443 are as in Fig. 1, in each 
case. Fluxes are given for an ISM propagation in a medium of n = 1 cm~ 3 , although 
again we have checked this is not a relevant parameter as discussed above. We find that 
clouds located from ~20 to ~30 pc produce an acceptable match to MAGIC data. In 
the case of a more impulsive accelerator, the VHE predicted spectra is slightly steeper 
than that produced in the continuous case at the same distance, so that it provides a 
correspondingly better fit to the MAGIC spectrum. Fig. 2 also shows, apart from MAGIC 
data, EGRET measurements of the neighborhood of IC 443. We recall that these two 
sources are not located at the same place, what we emphasize using different symbols. 
Fig. 2 shows that there is plenty of room for a cloud the size of that detected in front of 
IC 443 to generate the MAGIC source and not a co-spatial EGRET detection. In the case 
of Fermi, measurement of this region will allow us to constrain the separation between 
the SNR and the cloud, since for some distances a Fermi detection is also predicted. 
The existence of a VHE source without counterpart at lower energies is the result of 
diffusion of the high-energy CRs from the SNR shock, which is an energy dependent 
process leading to an increasing deficit of low energy protons the farther is the distance 
from the accelerator. 

To clarify our previous assertion, and since our solution to the diffusion-loss equation 
is a function of time, we show the evolution of the flux along the age of the SNR. In Fig. 
3 we show the integrated photon flux coming from the position of the giant cloud as a 
function of time above 100 MeV and 100 GeV in the impulsive case. Different qualities 
of the accelerator (impulsive or continuous) produce a rather comparable picture. At the 
age of the SNR (the time at which we observe) Fermi should see a source only for the 
closest separations. On the contrary, the integrated photon fluxes above 100 GeV present 
minimal deviations, and a MAGIC source is always expected. 

Fig. 2 also presents the results of our theoretical model focusing in the energy range of 
EGRET. There, the CR spectrum interacting with a local-to-the-SNR cloud is obtained 
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FIGURE 1. Current CR spectrum generated by IC 443 at two different distances, 10 (solid) and 30 
(dashed) pc, at the age of the SNR. Two types of accelerator are considered, one providing a continuous 
injection (black) and other providing a more impulsive injection of CRs (red). The horizontal line marks 
the CR spectrum near the Earth. The y-axis units have been chosen to emphasize the excess of CRs in the 
SNR environment. 



assuming an average distance of interaction of 3-4 pc. A few hundred M located at this 
distance (~700 M for the case of an impulsive, and ~300 M for a continuous case) 
produce an excellent match to the EGRET data, without generating a co- spatial MAGIC 
source. Concurrently with, e.g., Gaiser et al. (1999), we find that the lowest energy data 
points in the EGRET range are produced by bremsstrahlung of accelerated electrons, 
curves that for simplicity we do not show in Fig. 2. 

As spinoff of the constraints provided by the observed phenomenology (e.g., the 
molecular environment and the position of the y-ray sources) in the setting of this model, 
we find that D\q should be low, of the order of 10 cm 2 s _1 . By varying the diffusion 
coefficient and studying its influence in our results, we obtain that if the separation 
between the giant cloud and the SNR is >10 pc, an slower diffusion would not allow 
sufficient high energy particles to reach the target material; thus, the MAGIC source 
would not be there. On the other hand, if the separation between the main cloud and the 
SNR is < 10 pc, we would have detected an EGRET source at the position of the cloud, 
which is not the case. We then grasp the value of D\q at 1.5 kpc from Earth, combining 
MAGIC and EGRET observations. Such values of D\q are expected in dense regions of 
ISM such as the one we study (Ormes et al. 1988, Gabici & Aharonian 2007). 
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FIGURE 2. MAGIC and EGRET measurement of the neighborhood of IC 443 (stars and squares, 
respectively) as compared with model predictions. The top (bottom) panel shows the results for an 
impulsive (continuous) case. At the MAGIC energy range, the top panel curves show the predictions 
for a cloud of 8000 M located at 20 (1), 25 (2), and 30 (3) pc, whereas they correspond to 15 (1), 
20 (2), 25 (3), and 30 (4) pc in the bottom panel. At the EGRET energy range, the curve shows the 
prediction for a few hundred M located at 3-4 pc. The EGRET sensitivity curves (in red) are shown 
for the whole lifetime of the mission for the Galactic anti-centre (solid), which received the largest 
exposure time and has a lower level of diffuse y-ray emission, and for a typical position in the Inner 
Galaxy (dashed), more dominated by diffuse y-ray background. The Fermi sensitivity curves (in blue) 
(taken from http://www-glast.slac.stanford.edu/software/IS/glast_latperformance.html) show the 1-year 
sky-survey sensitivity for the Galactic North pole, again a position with low diffuse emission (solid), and 
for a typical position in the Inner Galaxy (dashed). 
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FIGURE 3. Integrated photon flux as a function of time above 100 MeV and 100 GeV, solid (dashed) 
lines correspond to the case of the cloud located at 10 (30) pc. The horizontal lines represent the values of 
integrated fluxes in the case that the CR spectrum interacting with the cloud is the one found near Earth. 
The vertical line stands for the SNR age. EGRET and Fermi integral sensitivity, consistent in value and 
color coding with those in Fig. 2, are shown. 

PREDICTIONS FOR FERMI 

The sensitivity of Fermi allows the observation of the previous scenario, but this obser- 
vation also depends on the spatial resolution of the different clouds. Fig. 4 shows the 
predicted energy spectrum for the two cloud systems when Fermi is not able to resolve 
them independently. The final spectrum is the result of the addition of the two spectra in 
Fig. 3. The figure also depicts the best fit of the model to a power law. Table 1. shows 
the different spectral index for the different energy bands. 

At high energies, we should see a morphological and a spectral change from the 
position of the cloud (i.e. the center of MAGIC J06 16+225) towards the center of IC 
443. At a morphological level, the lower the energy, the more coincident with the SNR 
the radiation will be detected. At a spectral level: sufficient statistics should show that 
the lower the y-ray energy the harder the spectrum is. Fermi observations may also be 
sensitive enough to detect the same cloud that shines at higher energy, which ultimately 
will allow to determine its separation from the remnant, if the diffusion coefficient is 
assumed -as we showed-, or viceversa. 

CONCLUDING REMARKS 

Here we have shown that MAGIC J06 16+225 is consistent with the interpretation of CR 
interactions with a giant molecular cloud lying in front of the remnant, producing no 
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FIGURE 4. MAGIC and EGRET measurement of the neighborhood of IC 443 (stars and squares, 
respectively) as compared with model predictions. The panel shows the results for an impulsive case 
where two curves have been added. This will be the scenario for those sources that Fermi can not resolve 
independently. At the MAGIC energy range, the curves show the predictions for a cloud of 8000 M Q 
located at 20, 25, and 30 pc. At the EGRET energy range, the curve shows the prediction for a few hundred 
M located at 3-4 pc. The orange region show the best fit of the model to a power-law spectrum, and such 
way how the spectral index can depend on the different energy bands. The Fermi sensitivity curves (in 
blue) (taken from http://www-glast.slac.stanford.edu/software/IS/glast_latperformance.html) show the 1- 
year sky-survey sensitivity for the Galactic North pole, again a position with low diffuse emission (solid), 
and for a typical position in the Inner Galaxy (dashed). 

TABLE 1. Mean values of the spectral indeces for the different energy ranges of 
the best fit to a power law of the model explaining the separation between MAGIC 
and EGRET sources. 

0.1<E<1 1<E<5 5<E<10 10 < E < 100 

Spectral index -1.91 ±0.01 -2.60 ± 0.03 -2.77 ±0.18 -2.66 ± 0.29 



counterpart at lower energies. We have also shown that the nearby EGRET source can 
be produced by the same accelerator, and that in this case, a co-spatial MAGIC source 
is not expected. In our model, the displacement between EGRET and MAGIC sources 
has a physical origin. It is generated by the different properties of the proton spectrum 
at different locations, in turn produced by the diffusion of CRs from the accelerator (IC 
443) to the target. Specific predictions for future observations can be made as a result of 
this model as we have shown for the Fermi case. 
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